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The exact solution for the two particle spectrum of the Heisenberg-Peierls one 
dimensional spin chain is given by working in the fermionic representation. The 
\ resulting equations for the eigenvalues are, in some sense, similar to those of the 

Richardson's solution of the BCS model and must be solved numerically. 
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Low dimensional spin systems have attracted a great deal of interest over the last years 
because, both, theoretical and experimental reasons. In fact, from the experimental point of 



view, there are actually many quasi one dimensional compounds such as the organic series 
(BCPTTF) 2 X with (X = AsF e , PF e ), the cuprate CuGe0 3 or the (VO) 2 P 2 7 compounds, 
as well as other series as the TTFMS^C^CF^)^ with M = Cu, Au, Pt, or Ni, which can be 
well described as spin-Peierls chains [lj. Furthermore, we have now the evidence that the two 
dimensional antiferromagnetism is one of the crucial components of the high temperature 
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superconductivity |2j, and we have also Heisenberg spin-ladder compounds which can be 
viewed as intermediaries between one and two dimensions, and also as between half-integer 
(S = 1/2) and integer (S = 0, 1) spin chains. 

From the theoretical side, these systems are interesting since their low dimensionality 
makes that their properties are strongly affected by quantum fluctuations and also because 
of its apparent simplicity. For example, low dimensional electronic materials are known to 
be very sensitive to structural distortions driven by the electron-phonon interactions which 
breaks the symmetry of the original (i.e. without electron-phonon interaction) ground state, 
resulting in a less symmetric but lower energy state, in which the electrons and ions are 
shifted, from their symmetric positions, in a regular manner that creates a periodic variation 
of the charge density called charge density wave. This is the famous Peierls instability which 
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opens a gap at the Fermi surface of the one- dimensional electronic chain, transforming what 
would be a metal (in absence of electron-phonon interactions) in a semiconductor. A similar 
effect is expected for one dimensional spin chains which are unstable against a dimerized 
spin-singlet ground state . 

One of the most widely used model for spin interactions is the generalized one dimensional 
Heisenberg Hamiltonian which can be written, in terms of the quantum spin operators 
Sj = (1/2) <Tj (where cr^ are the Pauli matrices at the site i), as 

M 

H = 4 J^+iStSj+i) (1) 

8=1 

and that, after some transformations, can be seen also as a one dimensional Hubbard model 
for electronic interactions. The case when the coupling Ji.i+i between nearest spins is taken 
as a constant is exactly soluble after Bethe and his celebrated Bethe ansatz [3] . Later it was 
solved again by applying the Quantum Inverse Scattering methods, with which this model 
results inside of a family of solvable models Q. However the Heisenberg model must be seen 
only as a first approximation to the real interactions, as long as a constant coupling between 
nearest-neighbors implies that the atoms are taken at fixed positions. If we account for the 
vibrations of the atoms, then the coupling between two nearest spins is proportional to the 
exchange integral and then it can be written, up to first order, as J n ,n+i = u — J\{v n — i> n +i) 
where u and J\ are constants and v n is the displacement of the nth atom from its equilibrium 
position. In the Born-Oppenheimer approximation, v n = (— l) n t>o (adiabatic limit) the 
Hamiltonian (JIJ can be written as 

M 

H = Y,(u + (-lTv)<T n cT n+1 , (2) 

71=1 

with v = 2J\Vq. The Pauli matrices er n act on the two-dimensional space rj n = C 2 whereas 
H acts on the 2 M -dimensional Hilbert space, 

n = ®n=lVm (3) 

which is the tensor product of M two-dimensional C 2 spaces. We use periodic boundary 
conditions on the lattice with M = even playing the role of the length of the chain. The 
total third component of the spin, 

M 

S3 = Y,rt (4) 

n=l 
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is a conserved magnitude that we shall use to describe the states of the system. In general, 
to find the exact spectrum of the Hamiltonian O) is a difficult task and several methods 

n n 

like bosonization p| and variational p have been applied and actually for v <C \u\, case 
of weak dimerization, the low energy region is believed to be equivalent to that of the 
Gaussian model pj. In this paper we want to find the low spin excitations over the state 
where all M spins are aligned and pointing in the — Z direction (note that the particular 
direction in which the spins are aligned is not relevant due to the rotational invariance 
of the Hamiltonian). To do this we shall fermionize the problem with the Jordan- Wigner 
transformation and use a method analogous to that of the Richardson's solution |s| of the 
BCS model of superconductivity. 



II. ONE PARTICLE STATES 



As said above, our starting point will be an eigenstate of H with all spins aligned in 
the — Z direction so that it is also an eigenstate of S3 with eigenvalue — M; starting from 
it, we want to find the eigenstates of the Hamiltonian that have one spin flipped in the 
+Z direction and we call them the one-particle states. The first task is to transform the 
Hamiltonian (J2J) in a fermionic Hamiltonian by means of the Jordan- Wigner transformation 
since, although this is not strictly necessary for these one-particle states, it will be need for 
the two-particle states, so we start by defining the fermionic operators: 

a t = K(l)af , (5) 
at = K{lW, (6) 

being of = (1/2) (erf ± iof) and 

K{l) = j^>t°J = l J\{-o*). (7) 

3=1 

These operators are true fermionic operators in the sense that |<4,a m j = 5^ m and 
|a|, a^j = = {ai,a m }. Furthermore, since 

of = 2a\ ai - 1, (8) 

the total third component of the states is, up to a constant factor, the number operator 
1 a\a,i, and we can identify the states with K spins flipped with those states that 
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have K particles of type a; furthermore, the total number of fermions is also a conserved 
quantity. Transforming from the a to the a operators, the Hamiltonian can be written as 

M / 1\ / 1 

H= 4 $>M) (M-g) { N ^~2 

M-l 

+ 2 ^ ci(u, v) (a\ai +1 + a[ +1 ai) - 2e lnN c M {u, v) (a^ai + a\a M ) , (9) 
i=i 

being iVj = aja/ and ci(u, v) = (u + (— Due to the last term, the Hamiltonian reads 
different depending on the number of particles that has the state on which H acts. Howewer, 
considering separately the subspaces of even and odd total number of fermions, we may write 
H in a simple universal way, that is 

M M 

H = uM + 4^ v) {NiN t+1 - N{) + 2 £ q(u, v) (a\a l+1 + a\ +l ai) (10) 

where 

a-M+i = —o-i, oL+i = ~ a i f° r states on which N is even 

a>M+i = a i, a L+i = a i f° r states on which N is odd, (11) 

then, acting on states that have an even number of fermions the Hamiltonian (|1U|) has 
antiperiodic boundary conditions, whereas the boundary conditions are periodic when H 
acts on states with an odd number of fermions. Trasforming to momentum space with the 
Fourier transformation given by 

i M 

bj = -^Ee^'o,, (12) 



1 M 

L _Y- e -ikW h (13) 



C'l 



where 



jj j if m is odd, 

ff (2j - 1) if m is even, j = 1, • • • , M 



AI 

and m means the number of fermions of the subspace on which H acts; then the operators 
fet and b n create or destroy, respectively, a fermion of momentum k(n) . With these trans- 
formations, the Hamiltonian (jlUj) in the momentum space is composed of four pieces, plus 
the identity times uM, 

H = H l + H 2 + H z + H A + uM (14) 
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where 



M 



H x = -AuJ2<n)b%, (15) 

n=\ 
M 

H 2 = -Aiv^sm(k(n))b ] n+M/2 b n , (16) 



n=l 

M 



and 



"3 ^ 2-~i C U n 1 +n 2 -n 3 u n: i u n 2 u n 1 , \± ' ) 

111,112,713 = 1 

V M 

Ha — —A — p^faO-fcCni^Lt lA L L /-ION 

ni,n2,n3=l 



e (n) = 1 -cos(k(n)). (19) 



The space of states will be generated now by a vacuum |0 > with zero excitations, over 
which act the operators b\ and hi that create and destroy excitations that we call particles 
with momentum k(l); those states are described by kets of the form l 2 , • ■ • , l m > that 
represents a state of m particles with momentums k(li), k(l 2 ), • ■ ■ k(l m ). Obviously, the 
total number of particles is conserved by every piece of the Hamiltonian and then we can 
calculated the energy of the eigenstates with a given number of particles. 

The spectrum of energies in the one particle subspace can be obtained easily due to the 
fact that Hi and H 2 are the only pieces that gives a non null contribution. Acting with H 
on the states \l > and \l + M/2 >, we have 

#i|Z> =-4ue(l)\l>, (20) 

M 

H 2 \l> = -Aw sin(k{l)) \l + — > , (21) 

=-4« £ (/ + f)| i + f >, (22) 

M M 
H 2 \l + — > =-Aivsin{k{l + —))\l>, (23) 

and H 3 \l> = H 4 \l> = H 3 \l + M/2 > = H A \l + M/2 > = 0. 

Then we see that the two dimensional subspace generated by the states \l > and \l + 4f > 
is invariant under the action of H and consecuently the eigenvalues of H 1 + H 2 are those of 
the matrix, 

\ 

(24) 



—Aue(l) Aiv sm(k(l)) 
-4ivsm(k(l)) -4ue(Z + f ) 
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and the total energies will be these eigenvalues plus uM, i.e., 



Ef = uM -4{u±^J(ucos(k(l)) 2 + (v sin(fc(Z)) 2 ), 1 = 1, 



M 



(25) 



where now k(l) = (21tc/M) as corresponds to the subspace with one particle. The corre- 
sponding eigenvectors are those of the matrix 



n. 



i(si ± u cos(fc(/))) 
y ±v sin(/c(/)) 

where s t = [u 2 cos 2 (fc(/)) + v 2 sin 2 {k{l))] 1 / 2 and nf = [2s 2 ±2s l ucos{k{l))] 1 / 2 . 



(26) 



III. TWO PARTICLES STATE CASE 



In this case, the pieces of the Hamiltonian H 3 and H 4 acting on the states are not zero 
and it does more difficult to find the solution. The states are described by the kets of the 
form | , Z 2 > with, < li < l 2 < M. We can observe that they are eigenstates of Hi 
and that the action of if 3 on them gives new states \l'i,l' 2 >, in both cases they fulfill the 
condition mod(li + l 2 , M) = mod(l[ + V 2 , M). In the same form, the parts H 2 and H4 mix 
the states l 2 > with states l 2 > being mod(h + l 2 , M) = mod(l[ + l 2 + M/2, M). In 
other words, Hi and conserve the total momentum whereas H 2 and H4 change the total 
momentum by an amount of 7r. Then, the subspace generated by states of the form 

|tf( m ) > = g(l,m-l)\l,m-l>, (27) 

i=Allpossibles 

\m(m + M/2) > = Y h(l,m-l + M/2)\l,m-l + M/2 >, (28) 

i=Allpossibles 

is invariant by the Hamiltonian, and we can reduce to this subspace the problem of finding 
the eigenvalues and eigenvectors of H. Note that in equations (j27M28J) the functions g(h, l 2 ) 
and h(li, l 2 ) are antisymmetric in their arguments, that is g(l\, l 2 ) = —g(h, h) and h(li, l 2 ) = 
—h(l 2 ,li), in the same way as they are the fermionic states \h,l 2 >= — \l 2 ,h >■ 



A. Case u = 



Let us begin with the case u = 0. Then the total Hamiltonian has only the alternating 
term so the equations are simpler and we can see beter the method. Working in units of v 
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(or v = 1) and taking 



E(m) = iij)(m) + ip(m + M/2), 
and demanding that (H 2 + H4) S(m) = EE(m) we obtain that 
£ [g(l,m- I) - g(m - Z,Z)] = 



(29) 



8 { 2B(m) sin [Jfe(Z - m/2)] + - sin [k(l)\ (h(l - M/2, m - I) - h(m -1,1 + M/2)) + 



+ sin [jfe(m - I)] (h(l, m-l + M/2) -h(m-l- M/2, /))} , (30) 



and 



E [h(l, m-l + M/2) - h{m - / + M/2, I)} = 8 {2A(m) cos [k(l - m/2)} + 

I sin [k(l)\ {g{m - I + M/2, 1 - M/2) - g(l - M/2, m-l + M/2)) + 



+ sin [k(m - I)} {g{l, m - I) - h(m - Z, I))} , (31) 



where A(m) and B(m) are defined as 



M m ) = j^J2 sin l k ( l ~ m / 2 )] 9{U m-l), 



B{m) = — Y,cos[k{l -m/2)] h(l,m - I + M/2). 



(32) 
(33) 



These equations (f37IJ) and (f3*Tj) can be solved for g and h. After a rather tedious calculus 
we have obtained the solutions for g(h, l 2 ), h(l\, l 2 ) and the eigenvalues E in the form 
... 8£M(Z! - Z 2 ) cos - Z 2 + M)/2)] 

h{lu h) ~ D^EMM) + 

UB{h - l 2 ) sin [k((h + l 2 + M)/2)} sin [k((h -l 2 + M)/2)) cos [k{{h -l 2 + M)/2)) 



g(h,k 



D 2 (E,h,l 2 ) 
8EB(k-l 2 )sm [k{{h-l 2 )/2)} 
Di(E, h, l 2 ) 

UA{h - l 2 ) cos [k({h + l 2 )/2)\ sin [k({h - l 2 )/2)\ cos [k({h - l 2 )/2)\ 



D 2 (E,h,l 2 ) 

where the functions Di(E, l±, l 2 ) are defined as 

Di(E,h,l 2 ) = E 2 -16\sm(k(h))-(-iysm(k(l 2 ))]\ (i=l,2) 



(34) 



(35) 



and the eigenvalues E depend only on the total momentum of the state m = l\ + 1 2 , and for 
each value of m satisfy the consistency equation 

sin 2 [k{l-m/2)\" 
, D 2 (E,l,m-l) , 



8E\ 2 cos 2 [k(l- m/2)]' 
~MJ I ^ 



L>i(£,Z,m-Z) 



E 



1. 



(36) 
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Notice that ()36|) is cuadratic in the energy E and hence it gives symmetric eigenvalues 
±E as corresponds to the charge conjugation symmetry of the Hamiltonian. Note that 
this consistency equation fffil) is of the type of the eigenvalue equation for the Richardson's 
solution for the BCS model [8j , so the study of the completeness of the solution is analogous 
to that of this model. 




-2 ; 
-3 ; 

--4 - 



Figure 1: Graphical solution of the equation (j.36j) for m = 3 and M = 12. Dashed line corresponds 
to the G(E) = (8E/M)~ 2 term and the continuus line is F(E). 

As an example, in Fig. ([TJ we plot the graphical solution of Eq. ()36|) for the case 
M = 12 and m = 3; there we plot the curves G(E) = (8E/M)~ 2 and F(E) = 
(j2i C °D}(mm~t) ) ^%^ETm^) * Here it is evident that the number of solutions is deter- 
mined by the discontinuities of the function F(E) along the positive axis. 

B. General case 

Now we return to the general case when u and v are ^ 0. The action of the different 
parts of the Hamiltonian on the states ([271- l2*8]l is given by, 

H x \il){m ) >= -iuY^{£{l)+£(m-l))g{l,m-l) \l,m-l >, (37) 
i 

H x \ij){m +M/2) >= -4M^(e(/) + e(m-/ + M/2)) 

i 
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h(l, m-l + M/2) \l,m-l + M/2 >, 

H 2 \ip(m ) >= -UvJ2(- sm(k(l)) g(l - M/2,m - I + M/2) 

i 

+ sin(fc(m - I)) g(l, m- l))\l,m - I + M/2 >, 

H 2 \i\){m +M/2) >=4™^(sin(£;(Z))/i(Z-M/2,m-Z) 

i 

+ sin(fc(m - I)) h(l, m-l + M/2)) \l,m-l >, 

H 3 \^j(m ) >= 16uA(m)^sm(k(l) - s{m)/2) |Z, m — I >, 

i 

H 3 \ip(m +M/2) >= 16uB(m) £ cos(/c(Z) - s(m)/2) \l,m-l + M/2) >, 

i 

H 4 \^(m ) >= l6ivA(m) £ cos(A;(Z) - s(m)/2) \l,m-l + M/2) >, 



H 4 \ip(m +M/2) >= lQivB(m) J2 sin(fc(Z) - s(m)/2) \l,m-l >, 

i 



being now 



A(m) 



M 



J2sm(k(l)-s(m)/2)g(l,m-l), 



B(m) = — J2cos(k(l) -s(m)/2)h(l,m-l + M/2), 



and 



s(m) 



2tx 
~M 



(m — 1). 



(38) 

(39) 

(40) 
(41) 

(42) 

(43) 

(44) 

(45) 
(46) 

(47) 



The sums are understood to all different values of the index and the integer labels of the 
states are module M. 

The energies will be obtained from the compatibility of the eigenvalue equation 

+ H 2 + H 3 + if 4 )(#(m) + *(m + M/2)) = E (*(m) + *(m + M/2)). (48) 

Equating the coefficients of every state in the first member and in the second member, we 
can observe that the relations between them mix only the g(l,m — l),g(l — M/2,m — l+M/2), 
h(l—M/2) and h(l—M/2, m—l) terms, besides of the functions A(m) and B(m) that depends 
of all coefficients. This relations are expressed as the system of linear equations, 



/ 



(El - M) 



\ 



g(l,m- I) 
g (l-f,m+f-l) 
h(l,m+f-l) 
\ h(l + f,m-l) J 



= 16 



1 sm(k(l) - s -^){uA{m) + ivB{m)) ^ 
- sm{k(l) - s -^)(uA(m) + ivB(m) 
cos(A;(Z) - s -&±)(-ivA(m) +uB(m)) 
-cos(Jfe(Z) - s -^)(-ivA(m)+uB(m)) 



(49) 
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where M is the matrix 



with 



M 



di(l,m) 

GZx(7-f,m) 



-ci(/,m) 
\^ -c 2 (l,m) 



c 2 {l,m) 
ci{l,m) 



Ci(/,m) 
-c 2 (l,m) 
d 2 (l } m) 




c 2 (Z,m) 
-ci(Z,m) 




d 2 (Z 



M 
2 ' 



/;/ 



(50) 



ci(/,m) 
c 2 (l,m) 
di(l, m) 
d 2 (l,m) 



Aivsm(k(m — I)) , 
Aiv sm(k(l)) , 
-4tt(e(/) + e(m - /)) , 

M 

-4u(e(l)+e(m -/ + —)). 



(51) 
(52) 
(53) 
(54) 



From the equations (J4U|) . we can get the solutions of g(Z, m — I) and /i(/, I — m + M/2) 
by making use of the inverse matrix of (EI — M) . The matrix elements of that matrix are 
function of the I, m, and the eigenvalues E of the Hamiltonian H = Hi + H 2 + H 3 + if 4 . If 
we call Si j to the matrix elements of (EI — M) _1 we obtain, 



g(l,m — l ) = 16(cos(/c(/) 



)(si,3 - si )4 )(-^A(m) + uB(m)) 



s[m 



+ sin(fc(Z) - ^- z )(si i i - si, 2 ))(«A(m) + iuB(m))) 



/i(Z, m — I + 



16(cos(fc(Z) 
s(m 



s(m) 



)(s3,3 - s 3>4 )(-wv4(m) + ttB(m) 



+ sin(fc(Z) - -V Z )(s 3 ,i - s 3)2 ))MM + ivB(m))) , 



where the s functions are 
1 



Sll 



1 

D 
1 



ci(Z,m) (i? — d 2 (l,m)) + 



c 2 (l,m) 



E-di(l 



M 



m) (E-d 2 (l,m))\ [E-d 2 (l 



M 



m, 



M 



«12 
Sl3 



—ci(l,m)c 2 (l,m) [2E - d 2 (l,m) - d 2 (Z — 2~' m > 

1 / 2 2 / M 

— ci(i,m) I ci(/,m) -c 2 (/,m) + I £ - g?i(Z - — ,m / 



E - d 2 (l 



M 

T' 



m 



1 / 2 2 / M \ 

S14 = —c 2 (l,m) (— Ci(/,m) +c 2 (/,m) + I £ - — 2~' m V {E - d 2 (l,m)) 
s 22 = ^(E-d 2 (l,m)) (c 2 (l,m) 2 + (E-d 1 (l,m)) - d 2 (Z - y , m))) + 



(55) 



(56) 
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1 ,2 ( ^ , n M \\ 

— ci(/,m) IE — d 2 (l — — , m) I 



1 / 2 2 ( . M 

■S23 = -^c 2 (/,m) (cl(/,m) -c 2 (l,m) - (E - d x {l , m)) IE — d 2 (l — — , m) 



524 : D 



— ci(i, m) [— ci(l,m) +c 2 (l,m) - (E — di(l,m)) (E — d 2 (l,m 



1 ,2 . 1 /„ . „ M 



S33 = -^ci(/,m) (J5-di(Z,m)) + — f J5-di(Z- — ,m' 



c 2 (/, m) 2 + - di(l, m)) ( i? - d 2 (l - 



M 
T 

s 34 = — ci(Z, m)c 2 (l,m) \ -2E + dx(l,m) + - — ,m)J 

1 , 2 /_, . M \ 

S44 = — ci(/,m) IE — cti(/ — — , m) I 

+ ±(E-d 1 (l,m)) (^c 2 (l,m) 2 +(E-d l (l-^,m)^ (E-d 2 (l,m))^j (57) 

being D the determinant of -EI — M, i.e., the characteristic polynomial of M 

D = ci(/,m) 2 (cl(/,m) 2 -c 2 (/,m) 2 + (E - d 1 {l,m)) (E - d 2 (l,m))) 

+c 2 (Z,m) 2 ^-ci(/,m) 2 + c 2 (/,m) 2 + - d x (l - —,m)j (E - d 2 (l,m))\ 

+ (c 2 {l, mf (E - di(l, m)) + (E-d 1 {l- m)) 

(ci(Z, mf + (E — di{l, m)) (E - d 2 {l, m)))) (e - d 2 {l - — , m)\ (58) 



With the last expressions, we can substitute g(l, m — I) and h(l, m — I + M/2) given by 
equations (|55j)- (|56j) in the definitions for A(m) and B(m) given in the equations (|45l4(j|) . If 
we call 

Clil (m, 25) = £ sin (*K0 - ^) oos(lf (0 - ^(« 1)3 - SM ), (59) 

Ci, 2 K £) = ^ £ sin(fc(0 - ^) 2 ( Slll - s lj2 ), (60) 

c 2 ,iK E) = ^E oos(*(I) - ^) 2 (s 3 ,3 - (61) 

= ^EM*(0"^) 008(^(0-^)^8,1-^), (62) 

Then we obtain the equations 

A(m) = ci y \{m, E)(—ivA(m) + uB (m)) + ci j2 (m, E)(uA(m) + ivB(m)) (63) 
B(m) = c 2j i(m, E)(—ivA{m) + uB{m)) + c 2)2 (m, E){uA{m) + ivB{m)). (64) 



12 



The compatibility of that system impose the condition, 



ivci 7 i(m, E) + wci^m, E) — 1 MCi 5 i(m, E 1 ) + ivc\^.{m, E) 
—ivc 2 ,i(m, E) + uc2 :2 (m, E) uc 2> x{m, E) + iuc 2i2 (m, - 1 



(65) 



this equation plays the role of the Bethe ansatz equations and it can be solved by algebraic 
or numerical methods for any 1 < m < M and an arbitrary number of sites M. 

The determinant D is the four degree characteristic polynomial of the matrix M that 
correspond to the Hamiltonian piece if 12 = H\ + if 2- It is factorized as the product 



where E{ are the eigenvalues of if 12. In some cases it can happen that some of this factors 
will be common to all numerators of the functions g and h involved in the same set of 
equations ()49|) . Then, the corresponding eigenvalue Ei of if 12 will be also eigenvalue of the 
total Hamiltonian, i.e., solution of (J63j) . Due to that, we must take care since that eigenvalue 
will not appear if we do the simplification in the expressions for the g and h functions. 

For small M the solutions of the equation can be compared with the solutions 
obtained after diagonalizing the total Hamiltonian in both, the site space or the momentum 
space. For example, for M = 6 we have that the dimension of the state space is 15 whereas 
the invariant subspace generate by the states with m and m + M/2 has 5 dimensions. The 
solutions can be obtained by algebraic methods in the two cases. Using the Mathematica 
application, for m = 1, we obtain in the two cases the solutions 



1 



D(E) = U( E - Ei) 



(66) 



i=l 



E 



E 







-8 m 



(67) 
(68) 




E 



(12)V3 



f{u,v) 



(69) 



(70) 



where f(u,v) is defined as 
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The solution E = —8u is also eigenvalue of Hi 2 . Numerical evaluations show that these 
solutions are real as correspond to an selfadjoint Hamiltonian. 

It is interesting to get the values of the energies for the special cases v — and u = 0. 
In the first case, we obtain obviously the spectrum of the XXX-Heisenberg model and the 
equations (f4*Hj) become very simple since Hi is diagonal and Hi is null. In the second case, 
u = 0, the equations (J44H) are simplified considerably too, since now the functions d are null. 
In both cases the solutions are, 

E v=0 = {0, -8 m, -4m, 2 (-5 + V5) u, -2 (5 + Vh) u}, (73) 
E u =o = {0, 0, 0, *V2v, -4V2v}, (74) 

that have been obtained from (|H7I - ITT]) by making use of the respective substitutions for u 
and v. 
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